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Abstract. Multi-dimensional constrained covariant density functional theories were developed recently. In 
these theories, all shape degrees of freedom fJ All deformations with even yi are allowed, e.g., P20, fin, [J30, 
1^40, P42, P44, and so on and the CDFT functional can be one of the following four forms: the meson exchange or 
point-coupling nucleon interactions combined with the non-linear or density-dependent couplings. In this contri- 
bution, some applications of these theories are presented. The potential energy surfaces of actinide nuclei in the 
(/?20,/322,/?3o) deformation space are investigated. It is found that besides the octupole deformation, the triaxiality 
also plays an important role upon the second fission barriers. The non-axial reflection-asymmetric fS i2 shape in 
some transfermium nuclei with A' = 150, namely 246 Cm, 248 Cf, 250 Fm, and 252 No are studied. 



o 

=5 



(N 

q 
o 



x 



1 Introduction 

"Shape" gives an intuitive understanding of spatial den- 
sity distributions of quantum many-body systems includ- 
ing atomic nuclei. For the description of the nuclear shape, 
it is convenient to use the following parametrization 
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where are the mass multipole operators. A schematic 
show of some typical shapes is given in Fig. Q] [1 1. The 
majority of observed nuclear shapes is of spheroidal form 
which is usually described by foo- Higher-order deforma- 
tions with A > 2 such as /^o also appear in certain atomic 
mass regions ||2). In addition, non-axial shapes in atomic 
nuclei, in particular, the nonaxial-quadrupole (triaxial) de- 
formation /?22 have been studied both experimentally and 
theoretically [3,4,5 1. There is no a priori reason to neglect 
the nonaxial-octupole deformations, especially the yfe de- 
formation |6,7j3). 

Furthermore, more shape degrees of freedom play im- 
portant roles in the study of potential energy surfaces of 
atomic nuclei. Particularly, various shape degrees of free- 
dom play important and different roles in the occurrence 
and in determining the heights of the inner and outer bar- 
riers in actinide nuclei (in these nuclei double-humped fis- 
sion barriers usually appear). For example, the inner fission 
barrier is usually lowered when the triaxial deformation is 
allowed, while for the outer barrier the reflection asym- 
metric (RA) shape is favored ll9lfT0][TT1[T2l[T3l. Nowadays, 
it becomes more and more desirable to have accurate pre- 
dictions of fission barriers also for superheavy nuclei [14. 
[T31IT6 17 18 19] . It is usually customary to consider only 
the triaxial and reflection symmetric (RS) shapes for the 
inner barrier and axially symmetric and RA shapes for the 
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Fig. 1. (Color online) A schematic show of some typical nuclear 
shapes: (a) sphere; (b) prolate spheroid; (c) oblate spheroid; (d) 
hexadecapole shape; (e) triaxial ellipsoid; (f) reflection symmet- 
ric octupole shape; (g) tetrahedron; (h) reflection asymmetric oc- 
tupole shape with very large quadrupole deformation. Taken from 
Ref. (TJ. 



outer one [15,20,21,22]. The non-axial octupole deforma- 
tions are considered in both the macroscopic-microscopic 
(MM) models [23] and the non-relativistic Hartree-Fock 
theories Il24l . 

In order to give a microscopic and self-consistent study 
of the potential energy surface with more shape degrees of 
freedom included, multi-dimensional constrained covari- 
ant density functional theories are developed recently 11251 
[261 . In these theories, all shape degrees of freedom /3a, de- 
formations with even /u are allowed, e.g., foo, P30, /3yi, 
Pao^PaIjPaa, and so on. In this contribution, we present two 
recent applications of these theories: the potential energy 
surfaces of actinide nuclei and the non-axial reflection- 
asymmetric /?32 shape in some transfermium nuclei. In Sec- 
tion[2j the formalism of our multi-dimensional constrained 
covariant density functional theories will be given briefly. 



EPJ Web of Conferences 



The results and discussions are presented in Section [3] Fi- 
nally we give a summary in Section |4] 



2 Formalism 

The details of the formalism for covariant density func- 
tional theories can be found in Refs. ll27ll28ll29l[30l[3T1[32l . 
The CDFT functional in our multi-dimensional constrained 
calculations can be one of the following four forms: the 
meson exchange or point-coupling nucleon interactions com 
bined with the non-linear or density-dependent couplings 
l26l . Here we show briefly the one corresponding to the 
non-linear point coupling (NL-PC) interactions. The Start- 
ing point of the relativistic NL-PC density functional is the 
following Lagrangian: 



L = <A (jy^ ~ M B ) Iff ~ -Clin - -£nl - -Cder - -£c 



(2) 
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are the linear, non-linear, and derivative couplings and the 
Coulomb part, respectively. Mb is the nucleon mass, as, 
ay, ars , <xtv< Ps, 7s ■> 7v< $s , 8v, Sts , and 6jv are cou- 
pling constants for different channels and e is the electric 
charge. p$, prs, jv, and jjy are the iso-scalar density, iso- 
vector density, iso-scalar current, and iso-vector current, 
respectively. The densities and currents are defined as: 
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PTS = 
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Starting from the above Lagrangian, using the Slater de- 
terminants as trial wave functions and neglecting the Fock 
term as well as the contributions to the densities and cur- 
rents from the negative energy levels, one can derive the 
equations of motion for the nucleons. Furthermore, for sys- 
tems with time reversal symmetry, only the time-like com- 
ponents of the vector currents ((6) and (|7} survive. The re- 
sulted equation for the nucleons reads 

/V/ = (or ■ p + p(M + S (r)) + V(r)) ifr, = e^u (8) 

where the potentials V(r) and S (r) are calculated as 

S = a s p s +Psp\ + 7sp\ + Ss&Ps 

+ (ctTSPTS + StS &PTS ) T3, (9) 
V - aypy + JvPy + SyApvW 

+ {ajyprv + 5tv&Ptv)t'}- (10) 



An axially deformed harmonic oscillator (ADHO) ba- 
sis is adopted for solving the Dirac equation [25 ,261 1331 . 
The ADHO basis are defined as the eigen solutions of the 
Schrodinger equation with an ADHO potential 03411351 . 
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where 



V B {z,p) = -M{o 2 p p 2 + a? z z 2 ), 
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is the axially deformed HO potential and a> z and co p are the 
oscillator frequencies along and perpendicular to z axis, 
respectively. The solution of Eq. dTTb reads 



a (rcr) = C a <p„Xz)K'(p)-^e im '* X sAcr), 
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where <p n Xz) and R™' p {p) are the HO wave functions, 
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Xs is a two component spinor and C a is a complex num- 
ber inserted for convenience. Oscillator lengths b z and b p 
are related to the frequencies by b z = 1 / -\JMa> z and b p = 

i/V^v 

These basis are also eigen functions of the z component 
of the angular momentum j z with eigen values K = mi+m s . 
For any basis state <P a (rcr), the time reversal state is de- 
fined as <t>a{rcr) = T<t> a (rcr), where T = i<jyK is the time 
reversal operator and K is the complex conjugation. Appar- 
ently we have K & - —K a and = n a . These basis form 
a complete set for expanding any two-component spinors. 
For a Dirac spinor with four components, 



ipi(ra-) 



Z„#<P„(«r) 



(15) 



where the sum runs over all the possible combination of the 
quantum numbers a = \n z ,n r ,mi,m s }, and f" and g" are 
the expansion coefficients. In practical calculations, one 
should truncate the basis in an effective way [25 261 1331 . 

We expand the potentials V(r) and S(r) and various 
densities in terms of the Fourier series, 



f(p,<p,z)= V fn(p,z)—=exp(ip.<p). 



(16) 



The nucleus is assumed to be symmetric under the V4 group, 
that is, for all the potentials and densities we have 



f(p, 7T + ip,z)= f(p, 2tt -<p,z) = f(p, <p, z). 



(17) 



Thus the components ^'s satisfy fy=fZ=fp and all the 
terms with odd yU vanish. The expansion Eq. (fT~6b can be 
simplified as 

f{p, ip, z) = fo(p, z) —== + V f n (p, z) — cos(2nip), (18) 
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are real functions of p and z. 

The total energy of a nucleus reads 

£totai = J" d 3 r jj] u 2 ^ (»■/»+ 0M) 
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where the center of mass correction E c m . can be calculated 
either phenomenologically or microscopically. 

The intrinsic multipole moments are calculated from 
the vector densities by 

Qa, = j d 3 rp v (r)r A Y Afl (Q), (21) 

where Y All {Q) is the spherical harmonics and r refers to the 
proton, neutron or the whole nucleus. 

The potential energy surface (PES) is obtained by the 
constrained self-consistent calculation, 



Nc 1 2 

E' = (H) + J]-C n ((Q n )-p„) 2 



(22) 



where H is the RMF Hamiltonian, Q n 's are the multipole 
operators to be constrained and N c is the dimension of the 
constraining space. 

Both the BCS approach and the Bogoliubov transfor- 
mation are implemented in our model to take into account 
the pairing effects. For the pairing force, we can use a delta 
force or a separable finite-range pairing force [36. 37]|38l. 
More details of the multi-dimensional constraint covariant 
density functional theories can be found in Refs. B25II261 . 



3 Results and discussions 

3.1 One-, two-, and three-dimensional potential 
energy surface of 240 Pu 

In Ref. Il25l . one- (1-d), two- (2-d), and three-dimensional 
(3-d) constrained calculations were performed for the ac- 
tinide nucleus 240 Pu. The parameter set PC-PK1 is used 1 39 
In Fig. |2] we show the 1-d potential energy curves (PEC) 
from an oblate shape with 020 about -0.2 to the fission 
configuration with P20 beyond 2.0 which are obtained from 




Fig. 2. (Color online) Potential energy curves of 240 Pu with var- 
ious self-consistent symmetries imposed. The solid black curve 
represents the calculated fission path with V4 symmetry imposed: 
the red dashed curve that with axial symmetry (AS) imposed, the 
green dotted curve that with reflection symmetry (RS) imposed, 
the violet dot-dashed line that with both symmetries (AS & RS) 
imposed. The empirical inner (outer) barrier height 6 cmp is de- 
noted by the grey square (circle). The energy is normalized with 
respect to the binding energy of the ground state. The parameter 
set used is PC-PK1. Taken from Ref. (25l . 



calculations with different self-consistent symmetries im- 
posed: the axial (AS) or triaxial (TS) symmetries combined 
with reflection symmetric (RS) or asymmetric cases. The 
importance of the triaxial deformation on the inner barrier 
and that of the octupole deformation on the outer barrier 
are clearly seen: The triaxial deformation reduces the inner 
barrier height by more than 2 MeV and results in a better 
agreement with the empirical datum; the RA shape is fa- 
vored beyond the fission isomer and lowers very much the 
outer fission barrier. Besides these features, it was found 
for the first time that the outer barrier is also considerably 
lowered by about 1 MeV when the triaxial deformation is 
allowed. In addition, a better reproduction of the empirical 
barrier height can be seen for the outer barrier. It has been 
stressed that this feature can only be found when the axial 
and reflection symmetries are simultaneously broken [25 1. 

In order to see how the PES of 240 Pu becomes unstable 
against the triaxial distortion, 2-d PES's from calculations 
without and with the triaxial deformation were compared 
in Fig. [3] [25 1 . When the triaxial deformation is allowed, 
the binding energy of 240 Pu assumes its lowest possible 
value at each (f$20,Pm) point. At some points, especially 
those around the two saddle points, non-axial solutions are 
favored than the axial ones. The inner barrier height is low- 
ered by about 2 MeV. About 1 MeV is gained for the bind- 
ing energy at the second saddle point due to the triaxiality. 
In the regions around the ground state and in the fission 
isomer valleys, only axially symmetric solutions are ob- 
tained. 

A full 3-d PES has been obtained for 240 Pu (25). In 
Fig.|4]are shown only five typical sections of the 3-d PES of 
240 Pu in the (J322,/3io) plane calculated around the ground 
state, the first saddle point, the fission isomer, the second 
saddle point, at a point beyond the outer barrier, respec- 
tively. The following conclusions were drawn by examin- 
ing these 3-d PES's [25]: (1) The ground state and the fis- 
sion isomer are both axially and reflection symmetric. The 
stiffness of the fission isomer is much larger than that of 
the ground state against both the P22 and fiio distortions. 
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Fig. 3. (Color online) Potential energy surfaces of 240 Pu in the 
(PiQtPia) plane from calculations (a) without and (b) with the tri- 
axial deformation included. The energy is normalized with re- 
spect to the binding energy of the ground state. The numbers in 
(b) show the values of fj 2 2 at these points. The fission path is rep- 
resented by a dash-dotted line. The ground state and fission iso- 
mer are denoted by full and open circles. The first and second 
saddle points are denoted by full and open triangles. The contour 
interval is 1 MeV. Taken from Ref. |25[. 
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Fig. 4. (Color online) Sections of the three-dimensional PES of 
240 Pu in the (fii2,Pio) plane calculated at B 2 q= 0.3 (around the 
ground state), 0.6 (around the first saddle point), 0.9 (around the 
fission isomer), 1.3 (around the second saddle point) and 1.6 (be- 
yond the outer barrier), respectively. The energy is normalized 
with respect to the binding energy of the ground state. The con- 
tour interval is 0.5 MeV. Local minima are denoted by crosses. 
Taken from Ref. [25 1. 
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Fig. 5. (Color online) The inner (B' f ) and outer (B°) barrier heights 
of even-even actinide nuclei. The axial (triaxial) results are de- 
noted by open (full) symbols. The empirical values are taken from 
Ref. 1 40 1 and represented by grey squares. Taken from Ref. [25 1. 



(2) The second saddle point which is close to /J20 = 1-3 ap- 
pears as both triaxial and reflection asymmetric shape. (3) 
The triaxial distortion appears only on the top of the fission 
barriers. 

From the investigation of the one-, two-, and three- 
dimensional PES of 240 Pu, we can learn a lot about the 
importance of different shape degrees of freedom in differ- 
ent regions of PES in actinide nuclei. These information 
could be useful in further systematic calculations. 



3.2 Inner and outer fission barriers of even-even 
actinide nuclei 

Guided by the features found in the 1-d, 2-d, and 3-d PES's 
of 240 Pu, the fission barrier heights were extracted for even- 
even actinide nuclei [25 1. The calculated values were com- 
pared with empirical ones recommended in RIPL-3 ||401 . 

As it has been shown previously, around the inner bar- 
rier an actinide nucleus assumes triaxial and reflection sym- 
metric shapes. Thus in order to obtain the inner fission 
barrier height we can safely make a one-dimensional con- 
straint calculation with the triaxial deformation allowed 
and the reflection symmetry imposed. In Fig.[5la) we show 
the calculated inner barrier heights B l { and compare them 
with the empirical values. It is seen that the triaxiality low- 
ers the inner barrier heights of these actinide nuclei by 
1 ~ 4 MeV as what has been shown in Ref. 1131 . In general 
the agreement of our calculation results with the empirical 
ones is very good with exceptions in the two thorium iso- 
topes and 238 U. Possible reasons for these disagreements 
were discussed in Ref. 11251 . 
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To obtain the outer fission barrier height B°, the situ- 
ation becomes more complicated because more shape de- 
grees of freedom play important roles around the outer fis- 
sion barrier. In Ref. 11251 . 2-d constraint calculations were 
made carefully around the second saddle points for even- 
even actinide nuclei. In the lower panel of Fig.[5]we show 
the results of outer barrier heights B° and compare them 
with empirical values. For most of the nuclei investigated 
here, the triaxiality lowers the outer barrier by 0.5 ~ 1 
MeV, accounting for about 10 ~ 20% of the barrier height. 
One finds that the calculation with the triaxiality agrees 
well with the empirical values and the only exception is 
248 Cm. From the calculation with the axial symmetry im- 
posed, the outer barrier height of 248 Cm is already smaller 
than the empirical value. The reason for this discrepancy 
may be related to that there are two possible fission paths 
beyond the first barrier ll25l . 

In Ref. Il25ll . it was also examined the parameter de- 
pendency of the influence of triaxiality on the outer fission 
barrier and the lowering effect of the triaxiality on the outer 
fission barrier was also observed when parameter sets other 
than PC-PK1 are used. 



3.3 Non-axial octupole shapes in N = 150 isotones 

Nowadays the study of nuclei with Z ~ 100 becomes more 
and more important because it not only reveals the struc- 
ture of these nuclei themselves but also provides significant 
information for superheavy nuclei [41.42. 43.44 1. One of 
the relevant and interesting topics is how to explain the 
low-lying 2~ states in some TV = 150 even-even nuclei. In 
these nuclei, the bandhead energy E(2~) of the lowest 2~ 
bands is very low 11451 . It is well accepted that the octupole 
correlation is responsible for it. For example, a quasiparti- 
cle phonon model with octupole correlations included was 
used to explain the excitation energy of the 2~ state of 
the isotones with N = 150 |46). In Ref. (47), Chen et 
al. proposed that the non-axial octupole i^-correlation re- 
sults in the experimentally observed low-energy 2~ bands 
in the N = 150 isotones and the reflection asymmetric shell 
model calculations reproduces well the experimental ob- 
servables of these 2~ bands. 

The non-axial reflection-asymmetric B32 shape in some 
transfermium nuclei with N - 150, namely 246 Cm, 248 Cf, 
250 Fm, and 252 No were investigated with the multidimen- 
sional constrained covariant density functional theory [48 1. 
The parameter set DD-PC1 is used 11491 . For the ground 
states of 248 Cf and 25() Fm, the non-axial octupole deforma- 
tion parameter B32 > 0.03 and the energy gain due to the 
Byi distortion is larger than 300 keV. In 246 Cm and 252 No, 
shallow /?32 minima are found. 

The triaxial octupole F32 effects stem from the cou- 
pling between pairs of single-particle orbits with Aj = Al = 
3 and AK = 2 where j and / are the total and orbit angular 
momenta of single particles respectively and K the projec- 
tion of j on the z axis. In Fig. |6] we show the proton and 
neutron single-particle levels near the Fermi surface for 
248 Cf as a function of B^ with B20 fixed at 0.3. It was shown 
that the spherical proton orbitals n2fi/2 and7rh'i3/2 are very 
close to each other [ 48 1 . This near degeneracy results in oc- 
tupole correlations. As seen in Fig. [6] the two proton levels, 
[521]3/2 originating from 2/7/2 and [633]7/2 originating 
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Fig. 6. (Color online) The single-particle levels near the Fermi 
surface for (a) protons and (b) neutrons of 248 Cf as a function of 
f} 32 with/? 2 o fixed at 0.3. 



from I/13/2, satisfying the Aj — Al = 3 and AK = 2 con- 
dition, are very close to each other at B20 = 0.3. Therefore 
the non-axial octupole T32 develops and with /S32 increas- 
ing from zero, an energy gap appears at Z = 98. Similarly, 
the spherical neutron orbitals v2gg/2 and vl 715/2 are very 
close to each other [48|. The neutron levels [734]9/2 origi- 
nating from 1 715/2 and [622]5/2 originating from 2^9/2 are 
also close to each other and they just lie above and be- 
low the Fermi surface. This leads to the development of a 
gap at TV = 150 with B^ increasing. The Y32 correlation in 
TV = 150 isotones is from both the proton and the neutron 
and for 248 Cf the correlation is the most pronounced [4-8 1 . 



4 Summary 

In this contribution we present some applications of the 
multi-dimensional constrained covariant density functional 
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theories in which all shape degrees of freedom defor- 
mations with even /i are allowed. The potential energy sur- 
faces of actinide nuclei in the (J320,/322,/33o) deformation 
space are investigated. It is found that besides the octupole 
deformation, the triaxiality also plays an important role 
upon the second fission barriers. For most of even-even 
actinide nuclei, the triaxiality lowers the outer barrier by 
0.5 ~ 1 MeV, accounting for about 10 ~ 20% of the barrier 
height. The non-axial reflection-asymmetric shape in 
some transfermium nuclei with N = 150, namely 246 Cm, 
248 Cf, 250 Fm, and 252 No are studied. Due to the interac- 
tion between a pair of neutron orbitals, [734]9/2 originat- 
ing from v/i5/2 and [622]5/2 originating from vgg/2, and 
that of a pair of proton orbitals, [521]3/2 originating from 
nfi/2 and [633]7/2 originating from ^'13/2, rather strong 
non-axial octupole F32 effects have been found for 248 Cf 
and 250 Fm which are both well deformed with large axial- 
quadrupole deformations, yS2o ~ 0.3. For 246 Cm and 252 No, 
a shallow minima develops along the P22 deformation de- 
gree of freedom. 
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